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(54) Procede de traitement d'image numerique pour I'extraction automatique d*objets en forme 
de rubans 

(57) Procede de traitement d'une image numerique 
(J Q ) pour ('extraction automatique d'objets en forme de 
rubans, incluant une operation de squeletisation com- 
prenant des etapes pour : 

former des images Iiss6es (Jj) a plusieurs echelles 
(aj), et dans chaque image tissue (J s ) respective- 
ment : 

extraire des bords des objets, extraire des pixels 
medians potentiels (fij P ) associes au centre (Qj) 
d'un cercle de rayon (Rj) proportionnel (k) a 
I'echelle (aj, tangent a des bords en un couple de 
pixels distincts (E v E 2 ), et associes a une mesure 
de conliance d'alignement du centre (H^ et des 
pixels du couple, et extraire des pixels medians Q M 
par une premiere selection des pixels medians 
potentiels de meme localisation ayant la mesure de 
conf iance maximale, et par une deuxieme selection 
des pixels medians potentiels restarrts ayant locale- 
ment une intensite maximale selon la direction d'ali- 
gnement 

et dans I'image numerique (Jmed) formee des 
pixels medians extraits, construire des squelettes 
des objets par un pistage. 
Application : Systeme d'imagerie medicale 



FI6.1 


Primed by Xerac (UK) Business Services 

2.16.2/3.4 


^ 1 EP0 8 

Description 

Linvention concerne un procede de traitement 
d'une image numerique pour I'extraction automatique 
d'objets substantiellertient en forme de rubans repre- 
sentes sur un fond. 

Llnvention concerne aussi un disposrtif mettarrt en 
oeuvre ce procede, combine a un systeme de formation 
•d'images medical es. 

Llnvention trouve particulierement son application 
dans I'industrie de fabrication des dispositifs a rayons X 
incluant des systemes de visualisation d'arteriogram- 
mes, y compris tes dispositifs dotes d'un mode de fonc- 
tionnement f luoroscopique. 

Un procede de traitement d'image pour la repre- 
sentation robust e d'objets est deja connu de la publica- 
tion intrtuiee "Robust Object Representation Through 
Object-Relevant Use of Scale" par Bryan MORSE, Ste- 
phen PI2ER et alii, dans "104, SPIE, vol.2167, Image 
Processing, 1994, pp.104-115". 

Cette publication d6crrt un procede de segmenta- 
tion d'image pour sSparer des objets et les representer 
dans un forme approprtee k leur suivi (identification, 
tracking). Le but du procede est de reussir une telle 
representation en s'affranchissant des deteriorations de 
Timage dues au bruit, au changement de grandisse- 
ment sans changement de resolution, au flou etc.. 

Cette publication definrt un premier concept appele 
"CORE" qui est le lieu des points situes au milieu de 
I'objet tel que mesure k des echelles proportionnelles a 
la largeur dudit objet 

Un point est drt srtue au milieu, ou encore sur une 
ligne mediane de I'objet, s'il repond a deux conditions. 
La premiere condition est qu'il faut qu'il existe au moins 
deux points de bord situes a une distance r (appelee 
demi-largeur) ou bien rayon de ce point median. La 
seconde condition est que la direction de ladite demi- 
largeur r doit §tre normale auxdits bords. 

Cette publication definrt un second concept pour la 
mise en oeuvre du procede, qui est que I'echelle de 
mesure, definie par un parametre a. doit etre propor- 
tionnelle k la demi-largeur r. 

Dans la publication crtee, le concept "CORE" ne 
represente pas des objets k une seule echelle d'un bout 
k I'autre d'une image, et m§me, ne represente pas un 
objet donne k une seule echelle a finterieur de I'objet 
lui-meme. Le concept "CORE" represente un objet 
simultanement dans une fourchette d'echelles formant 
des courbes separees k I'interieur d'un espace d'echel- 
les. 

Les etapes de ce procede de segmentation neces- 
siterrt d'abord le calcul des elements appeies "CORE* 1 
par des sous-etapes pour : 

1/ lisser I'image a plusieurs echelles pour construire 
un espace d'echelles qui decrit rimage a des multi- 
ples resolutions, 

2/ ^valuer a chaque echelle des nuages de points 
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de bords potentiels, qui torment un premier ensem- 
ble flou (fuzzy) appei6 "boundariness" et qui sort 
calculus au moyen d'operateurs appropries, 
3/ evaluer des nuages de points medians, qui for- 

5 ' merit' un second ensemble flou (fuzzy) appeie 
"medialness", qui sont definis par ('association de 
tous les points de bord potentiels. et qui sont calcu- 
lus selon les deux conditions decrrtes plus haut, 
dans I'espace muftiresolution, 

10 4/ trouver des aretes (ridges) dans I'ensemble flou 
appeie "medialness", qui sont des maxima d'inten- 
site de ce second ensemble flou. 

Le proced6 decrit dans le document cite impose de 
is local is er la, ou les lignes centrales d'un objet dans une 
image numerique en utiiisant des nuages de points de 
bords conduisant k des nuages de points medians pour 
garder la notion d'ensemble f lou le plus longtemps pos- 
sible, dans un espace m u I ti resolution, impliquant, pour 
20 chaque calcul, deux variables d'espace qui sont les 
coordonnees du point courant dans I'image, une varia- 
ble <j qui est le noyau des f iltres de lissage multiresolu- 
tion et qui est proportionneile k un rayon r a trouver, et 
une variable d'intensite liee audit point courant. Ces cai- 
25 oils doivent done etre menes a bien sur un nombre de 
points extrememerrt eieve, le tout dans une geometrie 
non Euclidienne. Ce precede n'est pas specialement 
applique k realiser la segmentation d'une image repre- 
sentant des objets d'une forme particuliere predetermi- 
ne? nee. 

Un but de la presente invention est de fournir un 
precede de traitement d'image numerique qui resout le 
probleme d' extraire, et non de segmenter, de maniere 
automatique, des objets substantial ement en forme de 
35 rubans representes sur un fond dans une image nume- 
rique, en mettant en oeuvre des etapes simples sur un 
petit nombre de pixels, menees a bien dans une geo- 
metrie Euclidienne, et peu consommatrices de temps 
de calcul. 

40 Ce but est atteint par un tel procede comprenant 
des etapes pour : 

1) former, k partir de ('image numerique, des ima- 
ges lissees a plusieurs echelles, et dans chaque 

45 image lissee respectivement : 

2) extra ire des pixels de bord des objets, 

3) extraire des pixels medians potentiels associes a 
une localisation d'un centre de cercle de rayon lie k 
I'echelle par une constante de proportionnaiite k, 

so tangent a des bords en un couple de pixels de bord 
distincts, et associes a une mesure de confiance 
que le centre du cercle et les pixels de bord du cou- 
ple sont substantiellement en alignement. 

4) extraire des pixels medians par une premiere 
55 selection des pixels medians potentiels extraits des 

differentes images lissees. qui orrt la mesure de 
confiance maximal e pour une m£me localisation, et 
par une deuxieme selection des pixels medians 
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potentiels restants qui ont localement une intensity 
maximale substantiellement selon la direction tfali- 
gnement, 

et dans llmage numerique formee des pixels 
medians extraits, * 5 
5) pister des pixels medians extraits pour construire 
des squelettes des objets. 

Le precede selon I'invention presente de nombreux 
avantages. Dans ce procede, les calculs sont effectues 10 
dans des images planes, traitees une par une. Ces cal- 
culs sont simples et le temps de calcul est court. Ce 
procede peut etre applique a rextraction et a I'identif ica- 
tion des vaisseaux dans des arteriogrammes comple- 
xes, a la caracterisation des stenoses ou des 15 
anevrismes dans des arteriogrammes, a I'eiimi nation 
d'objets ayant une geometrie connue egalement en 
forme de rubans de largeurs differentes de celles des 
vaisseaux, tels que les cdtes qui sont genantes dans I" 
examen d'arteriogrammes des poumons pour la detec- 20 
tion de I'emphyseme. 

Ce procede, et un dispositrf mettant en oeuvre le 
procede. sont decrits ci-apres dans des exemples de 
realisation illustres par les figures schematiques 
annexees, dont : 25 

la FIG.1A represente une image numerique origi- 
nale Jqo 

la FIG. 1 B est une image J 0 f iitree, 

les FIG. 1 C a 1 H sont des images Jfa) lissees a une 30 

meme echelle o\, et parmi ces Figs, 

la F1G.1C est une image numerique J iBG d'esti- 
mation du fond, 

la FIG.1D est une image Jj numerique dont le 35 
fond est extrait, 

la FIG. 1 E est une carte Kj binaire des pixels de 
bord, 

la FIG.1F est une image J^ed numerique des 
pixels medians, 40 
la FIG.1G est une image numerique J|so des 
pixels medians selectionnes par la selection 
isotrope, 

la FIG.1H est une image numerique J^-j- des 
pixels medians selectionnes par la selection 45 
anisotrope ou lateral e, 

la FIG.1 1 est une image numerique J SK des pixels 
medians chaTnes par pistage pour former le sque- 
lette des objets, 50 
la FIG.1 J est une image numerique reconstruite JR 
a I'echelle de ('image d'origine par integration des 
surfaces des cercles associes aux pixels medians, 
la FIG.2 represente les etapes du procede sous 
forme de blocs fonctionnels, 55 
les FIG.3A et 3B illustrent une etape d'extraction du 
fond a une echelle donnee aj et selon differentes 
orientations ; 


- la FIG.4 iilustre 1'etape d'extraction des bords a une 
echelle donnee a s ; 

- la FIG.5A iilustre dans I'etape de recherche des 
pixels medians potentieJs a une echelle donnee a i( 
'la generation d'un cerde de rayons Rj ; " 
la FIG.5B iilustre rextraction d'un premier pixel 
median potentiel a une localisation donnee associe 
a un premier couple de pixels de bord et la F1G.5C 
iilustre rextraction d'un second pixel median poten- 
tiel a la meme localisation, associe a un second 
couple de pixels de bord pour effectuer une sous- 
etape de selection isotrope ; 

la FIG.5D iilustre la sous-etape de selection laterale 
des pixels medians potentiels pour fournir les pixels 
medians ; 

la FIG.6 iilustre l'etape de pistage des pixels 
medians pour fournir les pixels du squelette ; 
la FIG.7 iilustre retape d'etude des variations des 
rayons d'un objet le long du squelette de cet objet ; 
la FIG.8 qui montre un dispositif medical a rayons X 
de formation d'images numeriques. 

PROCEDE 

Les figures sont decrites ci-apres en detail en rela- 
tion avec les differentes etapes d'un procede de traite- 
ment d'une image numerique, comprenant la 
representation d'objets en forme de rubans contrastant 
sur un fond. Ces etapes effectuent une squeletisation 
des objets de ('image, pour fournir de maniere automa- 
tique les pixels appartenant aux lignes medianes des 
objets et pour final ement reconstruire une image nume- 
rique ou les objets ont un niveau d'intensite contrastant 
plus vigoureusement vis-a-vis du fond, appele niveau 
d'intensite rehausse. 

Telle que representee sur la FIG.1 A dans ('applica- 
tion envisagee a titre d'exemple, I'image numerique 
d'origine J Q o es * un arteriogramme. Dans cet arterio- 
gramme, il existe des vaisseaux de diametres tres drffe- 
rents represerrtes par des rubans ici sombres sur un 
fond non uniforme. Les vaisseaux de fins diametres 
sont en particulier difficiies a discerner. 

Telle que representee sur la FIG.1J, apres le derou- 
lement du procede de segmentation, I'image J R est 
encore numerique et comprend exclusivement le sys- 
teme arteriel de ('arteriogramme de la figure 1 A, repre- 
sente avec un niveau d'intensite rehausse contrastant 
sur un fond uniforme. 

En reference avec la FIG.2, qui iilustre le procede 
sous forme de blocs fonctionnels, le procede de traite- 
ment d'image comprend de preference la suite d'etape 
decrite ci-apres : 

1/ Acquisition M00> d'une imaoe numerique Joo 

Telle que representee sur la FIG.1 A, I'image nume- 
rique est acquise par des moyens 100 de formation 
d'images numeriques, sous la forme d'une matrice bidi- 
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mensionnelle Jqo cte pixels ayant des niveaux d'inten- 
site numerises. Llmage acquise J 0 o represente, dans 
cet exemple. un arteriogramme et les vaisseaux tor- 
ment des rubans sombres sur un fond plus clair. 

.L'image acquise est alors soumise a un f iitrage 110 s 
qui rehausse I'intensite par exemple par une methode 
de tabulation (Look-up Table) classique, pour fournir 
une image numerique J Q representee a titre d'exemple 
sur la FIG.1B. 

w 

2/ Rltraqe multiresolution (10.1-10.N1 

L'image f iltree J Q est soumise a une serie de N fil- 
tres passe-bas 10.1 a 1 0.N ayant pour noyaux des fonc- 
tions gaussiennes isotropes d'ecart-types a-, k ct n is 
drfferents qui determined les resolutions correspon- 
dantes des N images J\ a J' N lissees obtenues par ce 
filtrage. Le nombre N est par exemple tel que 2 < N ^ 
100, et I'indice i est tel que 1 £ i ^ N. De preference N = 
10. Dans Pexemple ou N = 10, le filtrage multiresolution 20 
est effectue en passant sur l'image J Q , dix filtres gaus- 
siens 10.1 a 10.10 ayant respectvement des noyaux 
caracterises par des ecarts-type a 1 a a 10 , qui fournis- 
sent dix images numeriques J\ a J' 10 lissees a dix reso- 
lutions differ entes. Les ecarts-types <j\ sont designes 25 
par le terme "echelle af et les images correspondantes 
par J'j. 

Dans chacune des images de la serie, on etudie un 
type de vaisseau, et un seul, caract§rise par sa demi- 
largeur appelee Rj telle que : R s > a s . On etudie par 30 
exemple des vaisseaux de demi-largeur R , = k.a t ou k 
est une constante. De preference k = 3, d'ou Rj = 3 cr, . 
Le coefficient k et I'echelle a ( definissent done le type de 
vaisseau que Ton cherche a detecter dans chacune des 
images J'j de la serie. Le nombre et la valeur des echel- 35 
les CTj doivent etre choisis soigneusement afin de per- 
mettre la detection de tous les vaisseaux depuis les 
vaisseaux tres fins jusqu'au vaisseaux tres larges, ou 
autres objets en forme de rubans. 

40 

3/ Extraction du fond 

Une etape 20.1 a 20.N d' extraction du fond est tea- 
lisee respectivement sur chacune des images lissees J'j 
a I'echelle a,, pour eii miner les elements etrangers au 45 
systeme arteriel a chaque resolution. 

La FIG.3A represente un profil d'intenstte I le long 
d'une ligne L1 de pixels qui intersecte un vaisseau Va 
dans l'image lissee J'j. Ce profil d'intensite montre une 
partie Va1 a variations d'intensite lentes correspondant so 
au fond, et une partie Va2 formant un dome d'intensite 
correspondant au vaisseau. Uextraction du fond com- 
prend le passage sur ('image lissee J'j d'un f litre 20. i 
centre sur un pixel courarrt P(x,y) de la ligne L1 , et effec- 
tual une moyenne entre une intensite l A d'un pixel A ss 
srtue a une distance donnee D } a gauche du pixel cou- 
rarrt P(x,y) et une intensite l Q d'un pixel B situe a la 
meme distance D f a droite du pixel courarrt P(x«y) sur la 


ligne L1. Llntensite moyenne calcuiee est attribuee au 
pixel courarrt P(x,y). 

Par le passage de ce fiitre 20.i, le ddme d'intensite 
Va2 du au vaisseau est elimine, et une partie Va3 de 
niveau d'intensite est estim6e pour constrtuer le fond a 
Templacement du vaisseau. 

Cette etape de filtrage 20.i est effectuee en adop- 
tant, a I'echelle aj, une valeur de distance Dj telle que : 
Dj > Rj = k CTj ; par exemple Dj > 3 ct s . 

En reference a la FIG.3B qui represente schemati- 
quement des vaisseaux 2Va et 2Vb dans l'image lissee 
J'j, retape de filtrage 20.i est effectuee selon une plura- 
lite P de directions L1 a LP regulierement reparties 
dans le plan de l'image J'j de ti/P en k/P. De preference. 
P s 8 et on effectue le filtrage 20.i selon 8 directions 
espacees angulairement regulierement de ids en 7c/8. 

En chaque pixel courant P(x,y), la valeur d'intensite 
du fond est ensurte estimee comme la valeur la plus fai- 
ble de toutes les valeurs estimees en faisant varier 
I'orientation L1 a LP, a rechelle ctj donnee. Par cette 
operation de filtrage 20. i, une image estimee du fond 
Jj5G est obtenue pour chaque echelle o\, comme par 
exemple representee sur la FIG.1 C. 

Puis, comme par exemple representee sur la 
FIG. 1 D, une image Jj est formee en attribuant aux pixels 
des donnees d'intensites calcuiees en soustrayant, par 
une etape de soustraction 30.i, les vaieurs d'intensite 
estimees pour le fond et prelevees dans l'image J iBG . 
des donnees d'intensite de l'image lissee correspon- 
dante J'j. 

On obtient alors N images dont le fond a ete extrait, 
notees Jj, aux diverses echelles orj. 

Ce procede d'extraction du fond est tres avanta- 
geux parce qu'il ne cree pas de defauts, d'erreurs ou 
d'imprecision dans les regions de superposition des 
vaisseaux, et parce qu'il n'est pas gourmand en temps 
de cafcul. 

4/ Extraction des bords des obiets dans la serie 
d'imaqes J j 

Une etape 40. i est menee a bien pour extraire les 
bords des vaisseaux ayant des rayons Rj etudies dans 
chaque image Jj a une echelle donnee a it incluant une 
determination des gradients d'intensite au moyen du 
catcul des derivees d'ordre 1 de la fonction d'intensite 
dans l'image Jj, selon des axes de balayage perpendi- 
culaires notes x'x et y'y, en chaque pixel courant P(x,y) 
a traiter. 

De preference, le gradient d'intensite est obtenu 
par le passage, selon des directions parallel es' aux 
axes x'x et y'y. de filtres orientes connus sous le nom de 
"Sever et notes ci-apres Sovel x selon I'axe x'x et Sovel y 
selon I'axe y'y. La mise en oeuvre des filtres de Sovel 
est decrite dans un ouvrage d'enseignement general 
intitule "Computer Vision" par Dana H. Balard et Chris- 
topher M. Brown, publie chez PRENTICE-HALL, Inc.. 
Englewood, Cliffs, New Jersey, US. 07632 p. 77, t.6-M 
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en relation avec la FIG.3.10C. 

Comme illustr6 par la FIG.4, par ('application des fil- 
tres Sovel x et Sovel y on obtient des composantes 2 X et 
G y selon les axes x'x et y'y, du vecteur gradient dlnten- 
sjte G au pixel courant PJ^.y). Le calcul du rapport de s 4 
ces composantes G x et G y fournit la valeur de la tan- 
gente d'un angle 9 que fait ce vecteur gradient G par 
exemple avec Taxe x'x pris comme reference. Par le cal- 
* cul de la somme des projections des composantes S x 
et^Gy sur la direction 9, on deduit egalement le module 10 
flGfl, au signe pres t de ce vecteur gradient. On d§finit 
ainsi, en chaque pixel courant. le vecteur^ gradient par 
son angle d'orientation 9 et son module 

L'angle 9 est ensuite discretise avec un pas ti/M ou 
M est une constants Par exemple, la constante M est is 
choisie egale a 8 et le pas d'angle a 71/8. L'angle associe 
au module du gradient gGo est aiors non plus 9, mais un 
angle ©j qui est Tangle discret choisi le plus proche de 
9. Dans le cas d'un pas d'angle nJS, l'angle ©j associe 
au gradient prend les valeurs : 0 1 = 0 ; 0 2 = n/S ; © 3 = 20 
7t/4 ; © 4 = Zn/8 ; 0 5 = rc/2 ; 0 6 = SrdB ; ©7 = 3n/4 ; 0 8 = 
7*/8. 

Lextraction des bords dans une image Jj a une 
echelle donnee <jj comprend ensuite Termination des 
gradients non extremaux. Cette operation est realisee 25 
en balayant I'image traite J- t selon une plurality Q de 
directions regulieremerrt reparties dans le plan de 
I'image Jj, de n/Q en it/Q. 

Dans un exemple, on choisit Q = 8. Dans ce cas, les 
directions de balayage correspondent aux angles dis- 30 
crets ©! a © 8 . Les pixels courants montrant un gradient 
maximal par rapport a leurs voisins sur chaque direction 
donnee, sont s&ectionnes comme pixels de bord, et 
leur intensite est mise a zero, tandis que t'intensite des 
pixels montrant un gradient non maximal est mise a 1 , 35 
ou ('inverse. 

Dans un autre exemple, on choisit Q = 4. Par exem- 
ple, les directions de balayage correspondent aux 
angles discrets ©-,, © 3 , © 5 , ©7. Dans ce cas, lors du 
balayage, ies gradients ayant des angles ©2 et © 8 sont 40 
projetes sur la direction © 1( tandis que les gradients 
ayant des angles © 4 et © 6 sont projetes sur la direction 
© 5 . Les pixels courants montrant un gradient maximal 
sont selectionnes et extraits comme decrit plus haul 

D'a utres method es de balayage sont possibles 45 
selon un petit nombre d'axes regulierement angulaire- 
ment repartis dans le plan de (Image. Lors du balayage, 
les gradients sont projetes sur celle de ces directions la 
plus proche de son orientation calculee © 1 a © M et les 
pixels correspondants aux valeurs maximales de pro- so 
jections des gradients sont selectionnes comme pixels 
de bords. 

Une image binaire K it appetee carte des pixels de 
bords est construite a partir de ces donnees, a chaque 
echelle o\, comme par exemple morrtre sur la FIG.1 E. 55 


5- Extraction de pixels medians potentiels 
5.1 Localisation de pixels medians potentiels 

Cheque pixel de bord est muni de caracteristiques 
qui sont : 

ses coordonnees x,y dans I'image Kj 
et une direction ©j parmi a ©^j. 

La direction ©j est celle du vecteur gradient avec 
une imprecision irrferieure & kM. Dans I'exemple decrit, 
la direction ©j est celle du vecteur gradient avec une 
imprecision inferieure a nJS ; Tangle 9 du vecteur gra- 
dient est dans une fourchette de ± 71/8 centree sur la 
direction discrete ©j. 

Ces caracteristiques sont utilisees pour effectuer, 
dans chacune des images numeriques J i( dans une 
etape 50. i, une localisation de pixels medians poten- 
tiels, c'est-a-dire de pixels situes substantiellement a 
egal distance des bords des objets en forme de ruban, 
ou de vaisseaux. 

A cet effet, dans l'6tape 50.i, dans chaque image J; 
en reference & la FIG.5A, on gen^re d'abord un cercle, 
note Cj, ayant un centre Oj, un rayon de valeur Rj liee a 
T6chelle a s comme decrit precedemment, et forme de 
pixels Pj, P'j disposes a partir du centre Q a une dis- 
tance substantiellement egale a Rj, sur des rayons 
orientes de maniere discrete selon les directions discre- 
tes ©j des gradients d'intensite determtn§es plus haut. 
Puis, on balaye chaque image Jj de maniere conven- 
tionnelle, par exemple selon Taxe x'x, et, en reference k 
la FIG.5B, on applique en chaque pixel courant P(x,y) 
de localisation x,y, ie centre du cercle gen6re comme 
ci-dessus. On cherche s'il existe au moins un pixel El 
de bord substantiellement en coincidence avec un pixel 
Pj, P'j de cercle Cj sur un rayon dans une des directions 
discretes ©j, et si Indirection discrete d'orientation du 
gradient d'intensite G E1 assoctee a ce pixel de bord E-, 
coincide substantiellement avec la direction ©j de ce 
rayon. Un tel pixel de bord E 1 est alors situe sur un pre- 
mier rayon de tangence d'un bord incluant le pixel E 1 
vis-a-vis du cercle de centre localise en x,y et de 
rayon Rj. 

Si un tel premier pixel de bord est trouve, on 
cherche dans les m£mes conditions s'il existe un 
deuxieme pixel de bord, distinct du premier, substan- 
tiellement en coincidence avec un autre pixel du cercle 
Cj sur un rayon dans une autre des directions discretes 
©j et si la_direction discrete d'orientation du gradient 
d'intensite G K associee a ce pixel de bord E 2 coincide 
avec la direction de ce rayon. Un tel deuxieme pixel E 2 
de bord est alors situe sur un second rayon de tangence 
d'un bord B 2 incluant le pixel E 2 vis-a-vis du cercle de 
centre localise en x.y et de rayon Rj. 

Les pixels de bord E 1 et E 2 sont apparies si ces 
conditions sont reunies et si en outre les vecteurs gra- 
dients, soit pointent tous deux vers le centre Oj, soit en 
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sens oppose, selon que llmage J; a trailer est formee 
d'objets sombres sur fond clair ou I'inverse. Dans ces 
conditions, il existe un cercle Cj de rayon R| substantiel- 
lement tangent a un premier bord et a un second 
bord B2 respectivement.en Ej et E 2 . 5 

On peut ainsi d§finir dans la m£me image J;, ou 
dans differentes images Jj a differentes echelles aj, plu- 
sieurs couples de pixels E 1 , E 2 de bord correspondant a 
une localisation donnee x,y. Les centres Qj associes a 
ces drfferents couples sont appeles pixels medians 10 
potentiels f^p 

5.2 - Mesure de conf iance 

Dans E'etape 50. i, en outre, a chaque pixel median 15 
potentiel Qj P on associe une mesure de conf iance 
d'etre substantiellement dans I'alignement des pixels du 
couple E 1( E 2 . 

A cet effet. en reference aux FIG.5B et FIG.5C, on 
definit Tangle form6 entre les directions discretisees 20 
@ E1 , ®E2 des gradients d'intensite associes aux pixels 
du couple E 1t E 2 . et on definit une direction discrete pri- 
vilegiee y h , assoctee au pixel median potentiel Ojp qui 
est celle de la bissectrice de cet angle, et qui determine 
une mesure du demi -angle Ay entre les directions dis- 2s 
cretisees G^- L'indice h est une constante de dis- 
cretisation liee a j. La mesure de conf iance associee a 
un pixel median potentiel Ojp est donnee par le produrt 
de la plus petite des deux valeurs de gradient d'inten- 
srte, relatives respectivement aux pixels du couple E1 , 30 
E2, par la valeur du sinus du demi-angle Ay. Cette 
mesure de cortfiance est d'autant plus grande que le 
sinus du demi-angle Ay est proche de 1, done que le 
demi-angie Ay est proche de nJ2, correspondant au cas 
ou le pixel median potentiel est proche de I'alignement 35 
avec les pixels du couple E 1( E 2 , comme illustre par la 
FIG.5C. 

Les donnees relatives a la localisation des pixels 
medians potentiels Ojp et aux mesures de confiance 
associees sont alors extraites et stockees. Chacun des 
pixels 0\ P correspondants est un candidal pour former 
les pixels medians des objets. Les meilleurs candidats 
doivent etre selectionnes. Une image numerique J P a 
Techelle de depart de I'image J Q . formee des pixels 
medians potentiels Qp est illustree par la FIG.1F. 

6 - Extraction des pixels medians 

6.1 : Selection isotrope 

A une mime localisation donnee x,y, on examine 
dans une etape 60i s'il existe dans une meme image Jj 
ou dans differentes images J j( a differentes echelles, 
plusieurs pixels medians potentiels associes a des cou- 
ples de pixels de bord drfferents E 1( E 2 et associes a 
des mesures de confiance differentes. Dans ce cas, 
une premiere selection appelee isotrope est operee 
pour conserver le meilleur candidal. Cette selection 


consiste a determiner celui des pixels medians poten- 
tiels ayant cette localisation x,y dans une image Jj qui a 
en outre la mesure de confiance maximale, et a eliminer 
les autres pixels medians potentiels de meme localisa- 
tion dans la meme image Jj. Par aiileurs, on peut selec- 
tionner un pixel median potentiel ayant une mesure de 
confiance maximale vis-a-vis de pixels medians poten- 
tiels de m§me localisation x,y, mais extraits dlmages ^ 
a J M differentes. Dans ce cas, on favorise la detection 
d'objets ayant un diametre different de 2R it e'est-a-dire 
plus petit ou plus grand, ce qui est d'un grand interet en 
cas de croisement d'objets. Les pixels medians poten- 
tiels restants forment une image numerique J )SO telle 
que representee sur la F1G.1G. 

6.2 : Selection anisotrope 

Ces pixels medians potentiels restants, du fait de 
leur mesure de confiance maximale, sont substantielle- 
ment dans i'alignement des pixels du couple associe 

Les pixels medians potentiels restants onl des loca- 
lisations differentes. On effectue alors dans une etape 
70, une seconde selection, appelee anisotrope. A cet 
effet, I'image est balayee selon une droite perpendicu- 
laire a la direction privilegiee y h . Un pixel median poten- 
tiel Hjp est selectionne comme pixel median ftj M si son 
intensite est maximale vis-a-vis de ses voisins dans une 
fenetre inferieure ou egale au rayon R; defini plus haul. 

Une image numerique J MED telle qu'illustree par la 
FIG. 1 H est construite a I'echelle de I'image J Q avec tous 
les pixels selectionnes comme pixels mediants Oj M des 
objets en forme de ruban ou vaisseaux et qui ont ete 
extraits de toutes les images Jj. 

7. Pistage des pixels medians pour constituer le 
squelette des obiets 

Dans I'image J MEE > les pixels extraits dans I'etape 
precedente comme pixels medians des objets sont 
chatnes dans une etape 80 pour construire le squelette 
des objets, a partir d'un pixel median de depart, en 
determinant des pixels adjacertts. 

Les pixels medians extraits sont munis de teurs 
caracteristiques prelevees dans I'image J M ed et dans 
les images Jj d'ou ils ont ete extraits, qui sont : 

leur localisation x t y dans llmage J^ed 

leur intensite de depart dans I'image J MED 

la direction privilegiee y h qui leur est associee et 

qui a ete determined dans I'image lissee J}, 

I'echelle a\ de I'image lissee J s d'ou iis ont ete 

extraits. 

En reference a la FIG.6, (Image Jmed est balayee 
de facon conventionnelle. En un pixel median de depart, 
on cherche les meilleurs voisins parmi les autres pixels 
medians pour former le pixel predecesseur et le pixel 
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suiveur du squelette. 

A cet effet, on predetermine un angle ± by de 
recherche, par exemple ± it/8, autour de ia direction pri- 
vilegiee y h et on cherche les meilleurs voisins dans la 
zone .sectorielle* ainsi for'mee entre les droites et y n * s 
formant autour de vt/ h des angles ± 5y. Dans cette zone 
sectorielle, on predetermine une distance d maximaie 
de recherche, par exemple d = 15 pixels. On cherehe 
alors les pixels predecesseur et suiveur qui ont I'inten- 
site la plus proche de celle du pixel de depart dans la 10 
zone sectorielle et dans la fourchette de distance. En 
outre, on admet de chaTner non seulement les pixels 
medians qui ont ete extraits d'une meme image lissee a 
une echelle <jj donnee, mais aussi des pixels medians 
qui ont ete extraits d'autres images lissees aux echelles is 
les plus proches. On predetermine par exemple que Ton 
admet comme meilleurs voisins des pixels medians qui 
ont ete extraits dans les images lissees a une ou deux 
echelles superieures et inferieures. 

Les pixels medians trouves comme meilleur 20 
predecesseur et suiveur, aussi appeles adjacents, sont 
attaches au pixel median considered Puis de proche en 
proche, on cherche d'une part les meilleurs predeces- 
seurs et d'autre part les meilleurs suiveurs et on cons- 
truit ainsi, a partir de llmage numerique J MED , un 
chainage de pixels formant le squelette des objets 
constitue des pixels medians c'est-a-dire des pixels 
situes a egal distance des bords des objets. Lorsqu'on 
ne trouve plus de predecesseur ou suiveur pour une 
chaTne donnee, on arrete la chaTne, et on commence 
une autre chaTne a partir d'un autre pixel median 
comme pixel de depart. II en resulte que les squelettes 
des differents objets sont extraits independamment les 
uns des autres, comme il lustre par 11 mage J SK de la 
FIG. 11. 

On attribue par exemple, pour les applications ulte- 
rieures, des numeros q d'une suite continue a chaque 
pixel de chames de pixels medians. 

8. Applications 

8.1 Reconstruction d'objets dans l'image 

Dans cette application, on realise une operation 
supplemental 80 de reconstruction des objets en 
forme de ruban en integrant les surfaces des differents 
cercles de rayons R { lies aux differentes echelles a\, 
ayant pour centres les pixels du squelette, en attribuant 
a ces surfaces les intensites de ces centres. 

On obtient ainsi une image numerique J Rp comme 
illustre par la FIG.1J, des objets dont i'intensrte est 
rehaussee puisque Ton a selectionne comme pixels 
medians Q M et ulterieurement comme pixel du sque- 
lette, les pixels d'intensite extremale. Cette image peut 
etre soit visualisee, soit enregistree sur des systemes 
appropries. Un certain nombre de donnees concernant 
chaque pixel du squelette sont stockees. Ces donnees 
sont : 


son numero q, 

le rayon Rj du cerde correspondant a Techelle de 
llmage d'ou le pixel est extrait. 

8.2 Elimination d'objets genants dans llmage 

Par la methode decrite plus haut, les objets sont 
extraits indhnduellement puisque les pixels du squelette 
sont determines par la methode de pistage decrite qui 
elimine, comme ne faisant pas partie du meme objet, 
les pixels predecesseur ou suiveur dont la direction pri- 
vilegiee M/ h est tr °P eloigned de celle des pixels pistes. 
Ainsi deux objets qui se croisent ne peuvent etre detec- 
tes comme un seul objet : ils sont detectes comme des 
objets individuels. 

Cet avantage est utilise pour une operation 81 d'eli- 
mination des objets genants de Timage. Ainsi, il se pose 
un probleme dans les arteriogrammes thoraciques pour 
la recherche de remphyseme ou les c6tes sont des 
objets en forme de ruban qui genent I'examen des vais- 
seaux, egalement en forme de ruban, mais de diametre 
different. Les objets formes par les cotes sont elimines 
en extrayant les objets en forme de ruban du diametre 
correspondant aux cotes. 

8.3 Detection des stenoses ou des anevrismes 

Dans une autre application, des courbes de varia- 
tions des rayons Rj en fonction du numero "q" du pixel 
du squelette sont generees par une operation 82 qui 
permet de detecter facilement une variation abrupte du 
diametre de I'objet, comme montre FIG. 7. 

Sur de telles courbes, on detecte ainsi des retrecis- 
sements abruptes ST qui peuvent correspondre a des 
stenoses, ou des augmentations abruptes du diametre 
de I'objet qui peuvent correspondre a des anevrismes. 
On est raisonnablement sur qu'une augmentation de 
diametre est relative a un anevrisme sur un objet donne 
et n'est pas due a un croisement d'objets du fait que les 
objets sont extraits de maniere independante. 

Cela est different du resuitat d'un procede de seg- 
mentation selon lequel on ne peut pas separer deux 
objets qui se croisent ou qui se recouvrent. 

DISPOSITIF A RAYONS X 

En reference a la FIG. 8, a titre d'exemple, un sys- 
teme de radiographie numerisee comprend une source 
1 de rayons X ; une table mobile 2 pour recevoir un 
patient ; un dispositif intensrficateur d'image 3 couple a 
un tube video 4, qui fournit des donnees a un systeme 
5 de traitement numerique d'image comprenant un 
microprocesseur. Ce dernier comprend plusieurs sor- 
ties dont une sortie 6 couplee a un moniteur 7 pour la 
visualisation de l'image radiographique ou des images 
dlntensite. 

L'image radiographique numerisee peut contenir 
512 x 512 ou 1024 x 1024 pixels codes sur 8bits ou 10 
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bits. Chaque pixel peut etre ainsi affecte d'un parmi 256 
ou 1024 niveaux d'intensite. Par exemple, les regions 
sombres ont un niveau d'intensite faible, et les regions 
claires de 1'image ont un niveau d'intensite eleve. 

Llmage numerisee peut etra obtenue en mode 5 
f luoroscopique. L'invention peut trouver son application 
particulierement dans le trartement d'images arteriogra- 
phiques. 

, Le trartement de I'image numerique selon le pro- 
cede decrit plus haut est mis en oeuvre dans ses diffe- 10 
rentes etapes et operations dans le systeme 5. Les 
donnees sont stockees dans une zone memoire non 
representee et affichees si besoin sur le moniteur 7. 
Des moyens d'enregistrement non representes peuvent 
§tre en outre utilises. 15 

Revendications 

1 . Procede de trartement d'une image numerique pour 
('extraction automatique d'objets substantiellement 20 
en forme de rubans representes sur un fond, ce 
procede incluant une operation de squeletisation 
des objets comprenant des etapes pour : 

1) former, a partir de 1'image numerique (J Q ). 25 
des images lissees (JJ a plusieurs echelies 
(oj), et dans chaque image lissee (Jj) respecti- 
vement : 

2) extraire des pixels de bord des objets. 

3) extraire des pixels medians potentiels (Q P ) 30 
associes a une localisation d'un centre (Oj) de 
cercle de rayon (Rj) lie a I'echelle (q) par une 
constante de proportionnalite (k), tangent a 
des bords en un couple de pixels de bord dis- 
tincts (E 1t E2), et associes a une mesure de 35 
confiance que le centre {Q$ du cercle et les 
pixels de bord du couple sont substantielle- 
ment en alignement, 

4) extraire des pixels medians Q M par une pre- 
miere selection des pixels medians potentiels 40 
extraits des differentes images lissees, qui ont 

la mesure de confiance maximale pour une 
meme localisation, et par une deuxieme selec- 
tion des pixels medians potentiels restants qui 
ont localement une intensity maximale subs- 45 
tantiellemertt selon la direction d'alignement, et 
dans I'image numerique (Jmed) formee des 
pixels medians extraits, 

5) pister des pixels medians extraits pour cons- 
truire des squelettes des objets. so 

2. Procede selon la revendication 1, dans lequel 
I'etape 21 comprend Textraction d'un pixel comme 
pixel de bord d'un objet lorsque la valeur du gra- 
dient d'intensite associee a ce pixel est localement 55 
maximale dans une direction discretisee (@ v .. ©m) 
d'orientation de ce gradient, 

I'etape 3) comprend Textraction d'un pixel 


median potentiel C\ P associe a une localisation 
(x.y) du centre (fy du cercle sll existe un couple 
(E 1t E2) de pixels de bord distincts, substantielle- 
ment situes sur le cercie, sur des rayons de tan- 
.gence colncidant substantiellement avec des 
directions discretisees d'orientation des gradients 
des pixels de bord du couple, et associe a une 
mesure de confiance constituee par le produrt de la 
plus petite valeur des gradients d'intensite des 
pixels de bord du couple par la valeur du sinus du 
demi-angle entre les directions discretisees d'orien- 
tation des gradients correspondantes, et comprend 
I'extraction d'autant de pixels medians potentiels 
associes a une m§me localisation (x,y) qull existe 
de tels couples dans les differentes images lissees 
(Ji). 

3. Procede selon la revendication 2, dans lequel : 

I'etape 4) comprend la premiere selection 
d'un pixel median potentiel parmi plusieurs pixels 
medians potentiels de meme localisation dans les 
differentes images lissees (Jj) comme celui pour 
lequel le produrt constituant la mesure de confiance 
est maximal, et comprend la deuxieme selection 
d'un pixel median parmi des pixels medians poten- 
tiels restants comme celui qui presente une inten- 
site maximale consideree sur une direction 
d'orientation perpendiculaire a ia direction discrete 
(H*h) d'orientation de la bissectrice de Tangle entre 
les directions (©j) discretisees d'orientation des 
gradients des pixels de bord du couple (E 1( E 2 ) 
associe a ce pixel median. 

4. Procede selon la revendication 3, dans lequel : 

I'etape 5) comprend, pour le pistage des 
pixels medians extraits, I'association de chaque 
pixel median extrait avec des caracteristiques 
inciuant : 

sa localisation dans rimage numerique (Jmed). 
son intensite dans I'image numerique (Jmed). 
I'echelle (aj de I'image lissee (JJ d'ou il est 
extrait. 

une direction privilegiee (vh) constituee par la 
direction discrete d'orientation de la bissectrice 
de Tangle entre les directions (©j) discretisees 
d'orientation des gradients des deux pixels {E^ , 
E 2 ) de bord du couple qui lui est associe dans 
llmage lissee (JJ d'ou il est extrait, 
et la determination, a partir d'une localisation 
d'un pixel median extrait. d'un pixel median 
adjacent defint comme un pixel median extrait 
ayant des caracteristiques de distance (d), 
intensite (I), echelle (ai), et direction privilegiee 
(Vh), dans des fourchettes predeterminees, et 
le chatnage de proche en proche de pixels 
medians adjacents. 
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5. Procede selon Tune des revendications 2 a 4, dans 
lequel : 

I'etape 3) comprend la generation cfun cer- 
cle de centre (HJ localise au pixel courant P(x,y), et 
t forme de pixels disposes autour du centre k une 
distance substantiellement egale a la valeur de 
rayon (RJ proportionnelle a I'echelle (aj de II mage 
lissee (Jj), ces pixels du cercle etant en outre dispo- 
ses sur des rayons selon les directions discretisees 
(®i»©m) d'orientation associees aux gradients 
d'intensite des pixels de bord, 
et comprend le passage de ce cercle sur I'image lis- 
see, la fecherche d'un premier pixel de bord (E t ) 
coTncidant substantiellement avec un pixel du cer- 
cle, et ayant une direction discrete (0 E1 ) d'orienta- 
tion de gradient coincidant substantieilement avec 
le rayon correspondant du cercle, et si ce premier 
pixel (E1) de bord est trouve, la recherche d'un 
second pixel (E2) de bord coTncidant substantielle- 
ment avec un pixel du cercle distinct ayant une 
direction (©^ discrete d'orientation du gradient 
coincidant substantiellement avec le rayon corres- 
pondant du cercle, pour former un couple (E1, E2) 
de pixels de bord associe au pixel centre (Oj) for- 
mant un pixel median potentiel (Qj P ). 

6. Precede selon Tune des revendications 2 & 5, com- 
prenant dans I'etape 2) d'extraction des pixels de 
bord, la determination des orientations des gra- 
dients d'intensite, la discretisation de ces orienta- 
tions selon une pluralite de directions (01 k @M) 
angulairement reparties de nM en ic/M dans le plan 
de I'image lissee k I'echelle consideree (aj) et 
I'application d'un critere pour I' extraction des pixels 
de bord constitue par la condition qu'un pixel est 
detecte comme pixel de bord, et ulterieurement 
extrait, si la valeur du gradient est maximal e consi- 
deree local ement selon cette direction discrete 
determinee. 

7. Procede selon la revendication 6, comprenant dans 
I'etape 4) de seconde selection des pixels medians, 
la determination d'une direction d'orientation dis- 
cretisee de la bissectrice formant une direction pri- 
vilegiee (^h) associee au pixel median potentiel 
(Oj P ), le balayage de I'image lissee consideree per- 
pendiculairement k cette direction privilegiee, et la 
selection d'un pixel median d'intensite extremale 
selon cette direction dans une fenetre de largeur 
inferieure ou egale au rayon (Rj) du cercle. 

8. Procede selon Tune des revendications 1 k 7, com- 
prenant en outre une etape de reconstruction des 
objets en forme de ruban en integrant les surfaces 
de cercles de drfferents rayons (Rj) lies aux echel- 
les (aj) ayant pour centres les pixels du squelette, 
en attribuant k ces surfaces les intensites de ces 
centres. 


9. Procede selon Tune des revendications 1 a 8, com- 
prenant en outre une etape de determination des 
variations de diametres des objets en forme de 
ruban en examinant les variations du rayon (RJ en 

5 * fonction d'un numero (q) d'une suite attribue k cha- 
que pixel consecutif du squelette. 

10. Procede selon Tune des revendications 1 k 9, 
incluant une etape d'extraction du fond pratiquee 

10 sur les images lissees, avant I'etape d'extraction 
des pixels de bords, 
comprenant : 

une sous-etape de passage, sur les pixels de 
75 I'image lissee, d'un lot de P filtres reguliere- 

ment angulairement orientes de n/P en id?, qui 
chacun determine une moyenne d'intensite de 
pixels situes de part et d'autre du pixel courant, 
k une distance (Di) superieure au rayon (Rj) 
20 des vaisseaux dans ladite image lissee, et qui 

attribue cette moyenne au pixel courant, 
une sous-etape de formation cfune image du 
fond (Jibg) ou chaque pixel se voit attribuer la 
valeur d'intensite la plus basse parmi les inten- 
ds sites determine es par les filtres du lot, 

une sous-etape de soustraction du fond dans 
I'image lissee. 

11. Dispositif pour systeme medical, pour le traitement 
30 d'images medicales numeriques representant des 

objets en forme de rubans sur un fond, comprenant 
un systeme d'acquisrtion d'une image numerique, 
un systeme de traitement d'image ayant acces aux 
donnees de I'image numerique et k un systeme 
35 d'affichage des images acquises et les images trai- 
tees, dispositif dans lequel le systeme de traitement 
d'image comprend un microprocesseur pour mettre 
en oeuvre un procede selon Tune des revendica- 
tions 1 a 10. 
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FIG.3B 
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FIG. 5A 
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